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Abstract 



Singularities appear in numerous important mathematical models used 
in Physics. And in most of such cases singularities are involved in 
essentially nonlinear contexts. For more than four decades, general 
enough nonlinear theories of singularities have been developed. A 
critically important related feature is that, above certain levels in sin- 
gularities, the operation of multiplication, and in general, nonlinear 
operations on such singularities do inevitably branch in infinitely many 
ways, without the possibility for the existence of some unique natural 
or canonical way such nonlinear operations may be performed. Con- 
sequently, the choice in such branchings has to come from extraneous 
considerations. 



"We do not possess any method at all to derive 
systematically solutions that are free of 
singularities..." 

Albert Einstein : The Meaning of Relativity. 
Princeton Univ. Press, 1956, p. 165 
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0. Preliminaries 

Singularities have been present in Mathematics ever since the simple 
and natural looking issue of dividing by zero. And with the mathe- 
matisation of modern Physics, they are causing major difficulties in 
quite a number of disciplines of that field of science. Needless to 
say, various branches of engineering, as well as other fields of science 
and technology encounter similar difficulties due to singularities in the 
mathematical models employed. 

Here we shall consider singularities of functions / : X — > Y, where X 
and Y are some Euclidean domains or finite dimensional manifolds. 
Also, we shall consider singularities of various generalized functions 
which extend the above concept of function. In this way, the following 
two fundamental features of singularities will be of concern, namely 

• the extent of the set of singularities as a subset in the domain 
of definition of functions or generalized functions, 

• the behaviour of functions or generalized functions in the neigh- 
bourhood of singularities. 

The main conclusion obtained will be as follows : 

• In case no limitations are imposed on the above two features of 
singularities, as soon as multiplication, and in general, nonlinear 
operations are effectuated with generalized functions, there is an 
inevitable infinite branching in the way such operations can be 
defined. In other words, there is no canonical, natural or unique 
way multiplication, and in general, nonlinear operations can be 
defined for such singular generalized functions. Consequently, 
the specific choice of the result of multiplication, and in general, 
of nonlinear operations on generalized functions with singulari- 
ties has to be made based on extraneous considerations. 

A first systematic and far reaching mathematical approach to singu- 
larities was given in the 19th century by the theory of Functions of 
One Complex Variable. That was the time which led, among others 
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and not necessarily in a related manner, to the celebrated and yet un- 
solved Riemann Hypothesis, which shows the depth of the respective 
theory. 

As for singularities in the context of functions of one complex variable, 
one should recall the Great Picard Theorem, according to which an 
analytic function in the neighbourhood of an isolated singularity point 
that is an essential singularity will take on all possible complex values 
infinitely often, with at most a single exception. 

Consequently, in the neighbourhood even of one single and isolated 
singularity, one can expect a rather arbitrary behaviour, when one 
deals with more general functions than analytic ones. Not to mention 
the situation when the singularity points form a considerably larger 
subset in the domain of definition of a function. Therefore, the con- 
sideration of the above two aspects related to singularities is indeed 
appropriate, and in fact, necessary in the case of a deeper going and 
more wide ranging approach. 

Beyond the confines of analytic functions, it was the linear theory of 
Schwartz distributions which, starting with the mid 20th century, gave 
in certain respects a considerably more powerful and general treatment 
of singularities. 

The major limitation of the Schwartz approach is in its essential con- 
finement to linear situations, and its consequent inability to deal with 
singularities in a nonlinear context, and do so in a convenient and sys- 
tematic manner, without the recourse to what often are merely adhoc 
approaches. This contrast sharply with the remarkable natural ease 
and clarity the earlier complex function theory was applicable within 
a nonlinear context, restricting itself, as it did naturally, to analytic 
situations. 

Furthermore, as a rather unfortunate event, the more notorious than 
celebrated 1954 paper of Schwartz [1] claimed to prove that a nonlin- 
ear theory of distributions would altogether be undesirable, if not in 
fact, impossible ... 

Amusingly, this claim has attained a wider acceptance, [2], and con- 
sequently has for long distorted the perception of the situation con- 
cerning singularities within a general enough nonlinear context ... 
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As it happened, however, two general nonlinear approaches to singu- 
larities have nevertheless emerged, the first being the one listed under 
46F30 in the Mathematical Subject Classification of the American 
Mathematical Society, and with the second the Abstract Differential 
Geometry, or ADG, established by A Mallios in the 1990s, see [3-11] 
for some of the more relevant references in this regard. 

History, nevertheless, still seems to hang somewhat uneasily upon the 
issue of singularities ... 

Indeed, due to the remarkable generality, clarity and power of the 
Schwartz linear distributional approach, the issue of algebraic oper- 
ations on singularities has naturally and inevitably been restricted 
to addition alone, without a similarly general consideration of the 
operation of multiplication, an operation which ended up being in 
fact implicitly excommunicated form any possible suitable and gen- 
eral enough theory of singularities, in view of the mentioned 1954 
Schwartz paper and its long ongoing misinterpretations. 

And the surprising and hardly yet noted fact here is as follows. Al- 
though multiplication is closely related to addition, being in certain 
ways but a repeated addition, when it nevertheless comes to singular- 
ities, an essential difference appears between these two basic algebraic 
operations. Namely 

• addition can be extended from usual functions, be they regular 
or with singularities, to all sort of generalized functions, and such 
an extension appears to be naturally done in a unique, canonical 
manner, 

while on the other hand 

• due to most simple algebraic, more precisely, ring theoretic rea- 
sons, the extension of multiplication from usual functions, be 
they regular or with singularities, to generalized functions does 
no longer have such a naturally unique canonical way. 

Here is, therefore, the root of what is in fact no less than an infinite 
branching of the ways multiplications can naturally be defined for gen- 
eralized functions. 



4 



And no wonder, this root has so far been mostly missed due to the 
mentioned implicit omission to consider multiplication of singularities 
within a general enough, and not merely adhoc, context ... 

Within the general nonlinear treatment of singularities in 46F30, so 

far three classes of differential algebras of generalized functions have 

been used in a variety of problems, mainly for the solution of large 

classes of nonlinear systems of partial differential equations. 

It is instructive to recall the way these three classes of algebras relate 

to the above two fundamental features of singularities. 

In this regard, all these algebras are able to deal with the singularities 

the Schwartz linear theory of distributions can, since each of these 

algebras contains all the Schwartz distributions. 

The issue, therefore, is to what extent these algebras are able to deal 
with additional singularities. Let us then consider these algebras de- 
fined on any given Euclidean open set X C W 1 . Consequently, the 
generalized functions will extend various classes of usual functions 
/ : X — > E, where E is a finite dimensional real or complex vector 
space. 

Here it should be pointed out again that the class of admissible sin- 
gularities is large in no less than two significantly useful ways : 

• First, the singularities of the functions / : X — ► E considered 
can be given by arbitrary subsets E C X, subject to the only 
condition that their complementary X \ E, that is, the set of 
regular, or in other words, non-singular points, be dense in X. 
For instance, if X = R n is an Euclidean space, then the set E C 
X of singularities can be the set of all points with at least one 
irrational coordinate. Indeed, in this case the set X \ E of non- 
singular, or regular points is the set of points with all coordinates 
rational numbers, thus it is dense in X. A relevant and rather 
remarkable fact to note in this case is that the cardinal of the 
singularity set E is strictly larger than the cardinal of the set of 
non-singular points, namely, X \ E. 

• Second, there is no restriction on the behaviour of functions 
/ : X — > E in the neighbourhood of points in their singularity 
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Related to this second freedom in dealing with singularities, one should 
not forget its significant importance in applications. Indeed, as stated 
in the mentioned Picard Great Theorem, an analytic function in the 
neighbourhood of an isolated singularity point which is an essential 
singularity takes on all possible complex values infinitely often, with 
at most a single exception. 

Consequently, in the neighbourhood of a singularity, one can expect 
a rather arbitrary behaviour when one deals with more general func- 
tions than analytic ones. 

The first class of algebras of generalized functions was aimed to deal 
with singularities within a systematic and as widely applicable as pos- 
sible nonlinear theory, [12-29]. This class contains as particular cases 
all the subsequent classes of differential algebras of generalized func- 
tions constructed so far. Within this largest class, a special subclass 
- of so called nowhere dense algebras - was developed from the begin- 
ning, class which is able to deal with arbitrary closed nowhere dense 
subsets r C X of singularities, while no restrictions whatsoever are 
imposed on the behaviour of generalized functions in the neigbour- 
hood of singularities. 

The second class of algebras, [31-37], requires polynomial growth con- 
ditions on generalized functions in the neigbourhood of singularities. 
In this regard, these algebras of generalized functions - which are but 
a particular case of the infinite variety of all possible differential al- 
gebras of generalized functions introduced in [12-29] - suffer from a 
severe limitation. Namely, in the neighbourhood of singularities of 
their generalized functions, these algebras require a polynomial type 
growth condition, thus they cannot deal even with isolated singulari- 
ties such as essential singularities of analytic functions. 

The third class of algebras, [38-41, 9-11], is much more powerful than 
the class of so called nowhere dense ones, since these algebras are able 
to deal with arbitrary subsets £ C X of singularities, subject to the 
mild condition that the respective complementary subsets X \ E be 
dense in X, while again, no restrictions whatsoever are imposed on 
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the behaviour of generalized functions in the neigbourhood of singu- 
larities in E. 

An important fact to note here is that the subsets E of singularities 
can have a cardinal larger than that of the subsets X\E of nonsingular 
or regular points, since the condition that X \ E be dense in X can be 
satisfied even when X \ E is merely a dense countable subset of X, in 
which case E must of course be uncountable. 

As for the nonlinear operations on singularities, the nowhere dense 
algebras and those in the third class allow arbitrary smooth such op- 
erations, while in the second class only smooth operations with poly- 
nomial growth are possible. 

As a consequence of its restriction upon singularities, as well as upon 
operations on singularities, the second class of algebras cannot deal 
with a number of important problems which are easily treated within 
the nowhere dense algebras, or those in the third class. Among such 
problems are the following. 

The global version of the classical Cauchy-Kovalevskaia theorem for 
solutions of analytic systems of nonlinear partial differential equations 
cannot even be formulated, let alone solved, within the second class 
of algebras. 

On the other hand, the first class of algebras is already able to produce 
such a global version on the existence of solutions, [19,21,22]. 

Arbitrary Lie group actions, which are of major importance in the 
solution of partial differential equations cannot be defined within the 
second class of algebras. 

Here again, the nowhere dense algebras are already enough to define 
globally arbitrary Lie group actions. And as one of the consequences, 
one can for the first time obtain the complete solution of Hilbert's 
Fifth Problem, [28]. 

This problem, again, cannot be formulated, let alone solved, within 
the second class of algebras due to the polynomial mentioned type 
growth conditions which they require. 

Also, when defining differential algebras of generalized functions in 
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the case of domains X which are arbitrary finite dimensional smooth 
manifolds, the algebras in the first and third classes allow for consid- 
erably simpler constructions than those in the second class. 

At a deeper analysis, however, one that is done in terms of Sheaf The- 
ory, the essential difference between the nowhere dense algebras or 
those in the third class, and on the other hand, the algebras in the 
second class, is that the former are flabby sheaves, while the latter 
fail to be so, [9-11]. And as is known, [42], the lack of the flabbiness 
property in the case of spaces of functions or generalized functions is 
an essential indicator of their limitations in dealing with singularities. 

Lastly, it should be noted that, in [43], the study of a fourth class 
which is far larger then the above third class of algebras has been ini- 
tiated. 

1. Inclusion Diagrams and Reduced Power Algebras 
with the corresponding Ideals 

It is an elementary property of the linear vector space V(X) of Schwartz 
distributions that it can be represented as the quotient vector space 

(1.1) V'(X) = S 00 (X)/V 00 (X) 
of the vector subspaces 

(1.2) V°°(X) — ► S°°(X) — > (C°°(X)) N 

with the arrows " — >" representing usual inclusions " C" , and 

(1.3) S°°(X) = { s = (Vv)z,6N e (C°°(X)) N I s converges weakly mV'(X) } 

(1.4) V°°(X) = { v = {Xu)uen e S°°{X) | v converges weakly to mV'(X) } 

The remarkable fact, which has always been there in (1.1) - (1-4), is 
that in the right hand term of (1.2) we have the differential algebra 
(C°°(X)) N , yet in (1.2) one uses only vector subspaces of it. Indeed, 
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it takes little imagination to try to replace (1.2) with 

(1.5) 1(X) — > A(X) — >• (C°°(X)) N 

where -4.(X) is a subalgebra in (C°°(X)) N , while X(X) is an ideal in 
A(X), and thus instead of the quotient vector space in (1.1), obtain 
the quotient algebra 

(1.6) A(X) = A(X)/X(X) 

which may allow a nonlinear theory of generalized functions, thus of 
singularities as well. 

Indeed, for that purpose, it may be convenient to have the inclusion, 
that is, linear embedding 

(1.7) V{X) — > A(X) 

and of course, also suitable partial derivations on A(X), which in some 
convenient manner may extend the distributional partial derivations 
on V(X). Clearly, such partial derivations can easily be obtained on 
A(X), in case A(X) and I(X) are invariant under the natural term- 
wise partial derivations 

(1.8) (C°°(X)) N 3 s = (^„)„ 6N .— ► D?s = (DP^) um e (C°°(X)) N 
with p E N n . Namely, if one has 

(1.9) DP1(X)C1(X), DPA(X) CA(X), P eN n 

then one can simply define for p G N n , the corresponding partial 
derivation 



(1.10) A(X)3 s + X(X) i— >■ D p s + T(X) E A(X) 

Let us for the moment, however, deal only with the algebraic aspects 
of (1.5) - (1.7). An obvious immediate and simple way to obtain (1.5) 
- (1.7) would be to construct inclusion diagrams of the form, [12-29] 
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X{X) 



- (C°°PO) 



K 



(1.11) 



V°°(X) 



- S°°(X) 



which satisfy the condition 



(1.12) I(X)f\S°°(X) 



V°°(X) 



a condition which is both necessary and sufficient for the existence of 
the linear embedding 



(1.13) S 00 (X)/V 00 (X) 3 s + V°°(X) i— > s + X(X) G ^(X)/X(X) 



thus equivalently, for (1.7). 

Unfortunately however, inclusion diagrams (1.11) cannot be constructed 
in view of the simple fact that, [12-29] 

(1.14) ( V°°(X) . V°°(X) )f]S°°(X) g V°°(X) 

as simple counterexamples can show it. Indeed, it is easy to construct 
sequences v = {Xu)uen G V°°(X), such that v 2 G S°°(X), yet v 2 ^ 
V°°(X). For instance, when X — R, one can take Xv( x ) = cos(ux), 
and obtain indeed that v G V°°(X), v 2 G S°°(X), and furthermore v 2 
converges weakly to 1/2 in V(X), thus clearly v 2 V°°(X). 

Consequently, one can turn to the immediately more involved inclu- 
sion diagrams, [12-29] 
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X{X) 



A(X) 



(C°°(X)) r 



'1.15) V 



V°°(X) 



S°°(X) 



where V and S are vector subspaces, such that the following three 
conditions hold 

(1.16) I(X)f]S = V 

(1.17) V°°(X)f]S = V 

(1.18) V°°(X)+S = S°°(X) 

which, as it is easy to see, are both necessary and sufficient for the 
existence of the linear embeddings 

(1.19) V\X)S i— > s + V°°(X) G 5 00 (X)/V 00 (X) 

(1.20) 5/V9S + V i— >■ s + V°°(X) G 5 00 (X)/V 00 (X) 

(1.21) S/VBs + V i— )• s + I(I)GA(I)=i(I)/I(I) 

where the mappings (1.19), (1.20) are in fact vector space isomor- 
phisms. 

Now clearly, (1.19) - (1.21) give the desired linear embedding (1.7) 
of the Schwartz distributions into algebras of generalized functions, 
namely 
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(1.22) V(X) — >• A(X) = A(X)/X(X) 

In view of (1.5), (1.6), the algebras of generalized functions A(X) in 
(1.7), (1.22) are nothing else but reduced powers of the algebra C°°(X) 
of smooth functions on X. 

The general Model Theoretic, [44], construction of reduced powers, 
although hardly known as such among so called working mathemati- 
cians, happens nevertheless to appear in quite a number of important 
places in Mathematics at large. For a sample of them, one can note 
the following. The Cauchy-Bolzano construction of the field R of usual 
real numbers is in fact a reduced power of the rational numbers Q. 
More generally, the completion of any metric space is a reduced power 
of that space. Furthermore, this is but a particular case of the fact 
that the completion of any uniform topological space is a reduced 
power of that space. Also, in a rather different direction, the field *M 
of nonstandard real numbers can be obtained as a reduced power of 
the usual field R of real numbers. 

In view of the above, the use of reduced powers in the construction 
of algebras of generalized functions should not be seen as much else 
but a further application of that basic construction in Model Theory, 
this time to the study of large classes of singularities, and as such, to 
the solution of rather general nonlinear systems of partial differential 
equations, among others. 

As seen in the sequel, the ideals X(X) in (1.22) play the essential role 
in the inevitable infinite branching which occurs when multiplying 
generalized functions that have singularities above a certain level. 

2. Off-Diagonality Characterization 

A fundamental result in the nonlinear algebraic theory of generalized 
functions, see 46F30, is the simple and purely algebraic characteri- 
zation of the algebras of generalized functions (1.7) which are con- 
structed upon inclusion diagrams (1.15) - (1.21). In this regard, first 
we note that these inclusion diagrams can be simplified as follows. 
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In the inclusion diagrams (1.15) there are four spaces to be chosen, 
namely X(X),A(X), V and S. However, it is easy to see, [12-29], that 
such inclusion diagrams can be reduced to the simpler form depending 
only on two spaces (V, «S'), namely 



1(V,S') 



(2.1) 



V 



V°°(X) 



A(V,S') 



S°°(X) 



(C°°(X)) r 



U°°(X) 



where V,S' are vector subspaces in S°°(X), while 

(2.2) U°°(X) = {r H = {^,^,...)\^e C°°(X) } C S°°(X) 
is the diagonal in the Cartesian product (C°°(X)) N . 

As for the conditions (1-16) - (1.18), they now become 

(2.3) VfV = {0} 

(2.4) Z(V, SOPS' = W 

(2.5) V°°(X)®S' = S°°(X) 

The mentioned fundamental result regarding the characterization of 
algebras of generalized functions (1.7) constructed upon inclusion di- 
agrams (1.15) - (1.21) as simplified in (2.1) - (2.5) is the following : 



Theorem 2.1. [12-29] 



13 



Within a large class of inclusion diagrams (2.1) - (2.5), the conditions 
(2.3) - (2.5) are equivalent with the following off-diagonality one 

(2.6) l(V,S')f]U°°(X) = {0} 
Remark 2.1. 

It is both theoretically and practically important to note that, as seen 
in [12-29], there are infinitely many different inclusion diagrams (2.1) 
- (2-5). Moreover, they give infinitely many different corresponding 
algebras of generalized functions 

(2.7) A = A(V,S')/1(V,S') 

which, see (1.22), contain the vector space V\X) of Schwartz distri- 
butions. 



3. Inevitable Infinite Branching in the 
Multiplication of Singularities 

And now, we can come to the main issue in this paper, namely, to 
indicate the reason for the inevitable infinite possibilities in defining 
multiplication of generalized functions in case their singularities are 
above certain levels. 

The remarkable fact in this regard is that the respective reason is of 
a simple algebraic nature, namely, it is a direct consequence of the 
off-diagonality characterization in Theorem 2.1. above of the algebras 
of generalized functions (2.7 Constructed through the method of re- 
duced powers. 

Indeed, for a given pair V, S ' and a corresponding subalgebra A(V, S ') C 
(C°°(X)) N in an inclusion diagram (2.1) - (2-5), let us denote by 

(3.1) m(X,V,S',A(V,S')) 

the set of all ideals X in A(V, S ') which can occur in such inclusion di- 
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agrams (2.1) - (2-5). This means, therefore, that for every such ideal 
X G 1D(X,V,S',A(V,S')), there exists a corresponding algebra of 
generalized functions 

(3.2) A = A(V,S')/1 

which, in view of (1.22), contains the vector space V(X) of Schwartz 
distributions. 

Now, the essential point regarding the multiplication of singularities 
is that, in view of the inclusion 

(3.3) V®S' DU°°(X) 

in (2.1), it follows easily that the multiplication in each of the algebras 
of generalized functions A in (3.2) preserves the usual multiplication 
of smooth functions in C°°(X). 

On the other hand, regarding the multiplication of generalized func- 
tions that are not smooth - therefore, are elements in A\C°°(X) - it is 
well known that in general they do no longer preserve even the usual 
multiplication of continuous functions, this being one of the immedi- 
ate implications of the 1954 so called Schwartz impossibility result, 
[12-29]. 

Furthermore, in view of Remark 2.1. above, there are infinitely many 
ways according to which multiplication ends up being done, ways cor- 
responding to the various algebras A of generalized functions. And the 
possibility of this infinite branching of multiplication is manifested as 
soon as the generalized functions which are the factors in multiplica- 
tion belong to A \ C°°(X), and as such, and are farther and farther 
removed from C°°(X), or even form C°(X), that is, have a higher levels 
of singularity. 

Let us illustrate the above with a simple example. For that pur- 
pose, let us fix the pair V,S' and take as a corresponding subalgebra 
A(V,S') = (C°°(X)) N in the inclusion diagram (2.1) - (2-5). 
In such a case, one may expect a natural or canonical multiplication 
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if one would be able to single out in a suitable manner a certain ideal 
X G ID(X, V, S', (C°°(X)) N ), and thus obtain the corresponding alge- 
bra of generalized functions A in (3.2). 

Here however, the off-diagonality condition (2.6) interferes, leading to 
a rather involved structure for the set of ideals which satisfy that con- 
dition, as seen in 

Proposition 3.1. 

There is no largest ideal in the set 

(3.4) ID(X) 

of ideals X in (C°°(X)) N which satisfy the off-diagonality condition 

(3.5) lf]U°°(X) = {0} 
Proof. 

Let us again take X — R, together with v ' — (x ' u )uen, v " — (x eN £ 
(C°°(X)) N , where 

Xv( x ) — 1 + &m(i/x), x v( x ) — 1 + cos(vx), v<EN,xEX 
Then it follows easily that 

I' = v'(C°°(X)) N , 1" = v"(C°°(X)f e ID(X) 
However 

X = X' + X" 

is an ideal in (C°°(X)) N which fails to satisfy the off-diagonality con- 
dition (3.5), therefore 

(3.6) X i ID(X) 
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Indeed 

v = v > + v » e z 

and v = (Vv)j/eN, where 

tp v (x) = 2 + sin(ux) + cos(ux) > 0, v G N, x e X 
Consequently 

X = (C°°(X)) N 

thus it is not a proper ideal in (C°°(X)) N , and in particular, it does 
not satisfy condition (3.5). 

It follows that ID(X) does not contain ideals which may contain both 
ideals X' and X". 

□ 

In view of [12-29], it is obvious that there are infinitely many ideals 
in ID(X), and in fact, also in X e TD(X,V,S', (C°°(X)) N ). Further- 
more, various ideals in these sets clearly lead to significantly different 
multiplications in the corresponding algebras of generalized functions 
(3.2). 

4. Adhock Multiplications 

There is a considerable literature, [34], in which multiplications, and 
some more general nonlinear operations are performed upon a few spe- 
cific generalized functions with singularities, often including the Dirac 
delta distribution and possibly, its derivatives. Such studies, mostly 
undertaken by, and of interest to certain physicists and engineers, do 
not assume, let alone, aim to, or develop any general nonlinear theory 
for dealing with singularities, and instead, limit themselves clearly to 
a few special instances of such operations. 

The respective authors seem not to be cognisant of the wider existing 
relevant literature, let alone of the decades long existing systematic 
nonlinear theory of generalized functions. 
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5. The Apparent Novelty of the Inevitable Infinite 
Branching in the Multiplication of Singularities 

The remarkable survey of state of the art Physics, and also of a lot 
of the Mathematics underlying it in [45], see review [46], mentions in 
passing the Schwartz distributions. However, in its declared intention 
of guide to the laws of the universe, fails to point to the issue of non- 
linear operations on singularities, let alone, to mention the inevitable 
infinite branching which appears when multiplication or more general 
nonlinear operations on singularities are performed. 

On the other hand, most important physical phenomena and processes 
lead to singularities in their mathematical formulation. Furthermore, 
the respective mathematical models happen essentially to be nonlin- 
ear. One of them is the phenomenon of black holes in General Rela- 
tivity, a phenomenon to which the author of [45] has contributed in 
significant ways. 

And yet, when it comes to dealing with singularities in highly nonlin- 
ear contexts, there seems not to be a sufficient awareness about such 
critically important novel features as the inevitable infinite branchings 
which appear in such situations, when multiplication or more general 
nonlinear operations on singularities are performed. 
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